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Abstract 

A self-consistent analysis of K — > 2tt and K — > 3ir decays within a unique framework 
of chiral dynamics applied to the QCD-corrected weak nonleptonic quark lagrangian 
has been performed. The results on K —>■ 2n amplitudes at 0(p 6 ), including the 
value for e /e, are compared with the experiment to fix phenomenological B-factors 
for mesonic matrix elements of nonpenguin and penguin four-quark operators. The 
dependence of B-factors on different theoretical uncertainties and experimental er- 
rors of various input parameters is investigated. Finally, we present our estimates 
at 0(p 6 ) for the CP-asymmetry of linear slope parameters in the K — > 3n Dalitz 
plot. 
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The starting point for most calculations of nonleptonic kaon decays is an effective weak 
lagrangian of the form [T|. 

Cl(\AS\ = 1) = V2G F V ud V: s J2QO l} (1) 

i 

which can be derived with the help of the Wilson operator product expansion (OPE) from 
elementary quark processes, with additional gluon exchanges. In the framework of pertur- 
bative QCD the coefficients Ci are to be understood as scale and renormalization scheme 
dependent functions. 

In ([[]), Oi are the four-quark operators, defined either by combinations of products of 
quark currents (i = 1,2,3,4, non-penguin diagrams) or, in case of gluonic (i = 5,6) and 
electro-weak (i = 7, 8) penguin operators, by products of quark densities. The operators 
Oi with i = 1,2,3,5,6 describe weak transitions with isospin change AI = 1/2 while 
the operator O4 corresponds to AI = 3/2 transition and operators C^s - to mixture of 
AI = 1/2 and A J = 3/2 amplitudes. 

The general scheme of the meson matrix element calculation by using of the weak 
lagrangian (0) is based on the quark bosonization approach [Q. The bosonization procedure 
establishes the correspondence between quark and meson currents (densities) and product 
of currents (densities). Finally, it leads to the effective lagrangian for nonleptonic kaon 
decays in terms of bosonized (meson) currents and densities: 

nn/ T/v 5 U a fl ^ ja(mes) 1 5\ \ a , ja(mes) 
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The meson currents/densities J1/ R ^ and J^/r are obtained from the quark determinant 
by variation over additional external sources associated with the corresponding quark cur- 
rents and densities ||. From the momentum expansion of the quark determinant to 0(p 2n ) 
one can derive the strong lagrangian for mesons C e ff of the same order and the correspond- 
ing currents and densities J1/ Rll and J^m to the order 0{p 2n ~ l ) and 0(p 2n ~ 2 ), respectively. 
Thus, the bosonization approach gives us the correspondence between power counting for 
the momentum expansion of the effective chiral lagrangian of strong meson interactions, 

r(mes) _ r{p 2 ) , ^(p 4 ) , r{p 6 ) , 

*-s *-'.s ~ A -'.s ~ '-'s ' ■ ■ ■ 1 

and power counting for the meson currents and densities: 

£<f n ) => jjp"- 1 ) (currents) ; Cf ] J (p "" 2) (densities) . 

Some interesting observations on the difference of the momentum behavior of penguin 
and non-penguin operators can be drawn from power-counting arguments. The leading 
contributions to the vector currents and scalar densities are of Oip 1 ) and 0(p°), respectively. 
Since in our approach the non-penguin operators are constructed out of the products of 
currents J£ M , while the penguin operators are products of densities J£, the lowest-order 
contributions of non-penguin and penguin operators are of 0(p 2 ) and O(p ), respectively. 
However, due to the well-known cancelation of the contribution of the gluonic penguin 
operator 5 at the lowest order J4j, the leading gluonic penguin as well as non-penguin 
contributions start from 0(p 2 ) Q Consequently, in order to derive the currents which 
contribute to the non-penguin transition operators at the leading order, it is sufficient to 
use the terms of the quark determinant to 0(p 2 ) only. At the same time the terms of 
the quark determinant to 0(p 4 ) have to be kept for calculating the penguin contribution 
at 0(p 2 ), since it arises from the combination of densities, which are of 0(p°) and 0(p 2 ), 
respectively. In this subtle way the difference in momentum behavior is revealed between 
matrix elements for these two types of weak transition operators; it manifests itself more 
drastically in higher-order lagrangians and currents. 

This fact makes penguins especially sensitive to higher order effects. In particular, the 
difference in the momentum power counting behavior between penguin and non-penguin 
contributions to the isotopic amplitudes of K — > 3n decays, which appears in higher orders 
of chiral theory, leads to the dynamical enhancement of the charge asymmetry of the Dalitz- 
plot linear slope parameter |5], |3|] . 

In our approach the Wilson coefficients Cj in the effective weak lagrangian ([!]) are treated 
as the phenomenological parameters which should be fixed from the experiment. They are 
related with the Wilson coefficients Cf D (/•«), calculated in perturbative QCD ||, via the 
-Bi-factors: 

The coefficients Cf CD ([i) contain a small imaginary parts and can be presented in the form 
of a sum of z and y components 

C? CD ^)=Cr\») + rC\v\^ r = -|4|. (2) 

V udV us 

2 There is no cancellation of the contribution of the electromagnetic penguin operator Og at the lowest 
order and the leading contributions start in this case from 0(p ) 
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Respectively, the amplitudes of nonleptonic kaon decays also contain z- and y- components, 



A = A {z) + rA {y) , (3) 

which can be expressed in terms of z- and y-components of Wilson coefficients. The dom- 
inating contributions A® of the four-quark operators Oj and 5j-factors may be written 

as 

A^ = [ - Cl z ' y) (fi) + Ci z ' y) (fi) + C^ y) (/*)] B x (/,) JF> 

+Ci z ' y) (fi)B^) A® + Ct y \v)BM ^ (5) + C^\fi)BM A® . 

The observable effects of direct CP-violation in the nonleptonic kaon decays are caused 
by the ^-components of their amplitudes ([|). In particular, the ratio e fe can be expressed 
as 

/£\ UJ A^ 

Re(-)=lmX t (P -P 2 ), P I= 

\e/ V2\e\\V ud \\V us \ A)' 

where ImAj = ImV t *V tl i = \ V u b\\V c b\sm5; uo = A^/Aq , and Aj are isotopic amplitudes of 
K — > 2tc decay corresponding to the tttt final states with the isospin 7 = 0,2. 

There are the following theoretical uncertainties which appear both from short distance 
(Wilson coefficients) and long distance (effective chiral lagrangians and ^-factors) contri- 
butions to the kaon decay amplitude: 

• Dependence on the parameter Im r ~ Im X t arising from the imaginary part of Wilson 
coefficient (0). 

• Regularization scheme dependence which arises when calculating Wilson coefficients 
beyond the leading order (LO) of QCD in the next-to-leading orders: naive dimen- 
sional regularization (NDR), t'Hooft-Veltman regularization (HV). 



Dependence of Wilson coefficients on choice of the renormalization point a and QCD 

MS 



scale A^L = (325 ±110) MeV. 



• Dependence of meson matrix elements on the structure constants L 2 , L 3 , L 4 , L 5 , L 8 
of the general form of the effective chiral lagrangian introduced at 0(p 4 ) by Gasser 
and Leutwyler 0. 

• Dependence of meson matrix elements on the structure constants of the effective chiral 
lagrangian at 0{p 6 ) ||. 

• Factors Bi (i = 1,4,5,8) for dominating contributions of four-quark operators Oi to 
the meson matrix element (£|). 

• Dependence on choice of the regularization scheme to fix the UV divergences resulting 
from meson loops. 

In the present paper we have combined a new systematic calculation of mesonic matrix 
elements for nonleptonic kaon decays from the effective chiral lagrangian approach with 
Wilson coefficients C® CD ([i), derived by the Munich group for /x = 1 GeV and m t = 170 
GeV. For the parameter lm\ t we have used the result obtained in |J: ImA t = (1.33 ± 
0.14) ■ 10~ 4 . We performed a complete calculation of K —>■ 2n and K — > 3tt amplitudes at 
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0(p 6 ) including the tree level, one- and two-loop diagrams. For the structure coefficients 
Li of the effective chiral lagrangian at 0(p 4 ) we used the values fixed in [10 from the 



phenomenological analysis of low-energy meson processes. The structure coefficients of the 
effective chiral lagrangian at 0(p 6 ) have been fixed theoretically from the modulus of the 
logarithm of the quark determinant of the NJL-type model (see H for more details). The 
superpropagator regularization has been applied to fix UV divergences in meson loops. The 
isotopic symmetry breaking (tt -^--)] mixing) was taken into account at the tree level. In 
[ |TTf one can find more technical details of the calculation of K — > 2n amplitudes. 

In our phenomenological analysis the results on the K — > 2n amplitudes, including the 
value for e /e, are compared with the experiment to fix the phenomenological Sj-factors 
for the mesonic matrix elements of nonpenguin and penguin four-quark operators. As 
experimental input we used the experimental values of the isotopic amplitudes Aq2 fixed 
from the widths of K — > 2ir decays and the world average value Ree /e = (19.3 ±3.8) x 1(T 4 
which includes both old results of NA31 |T2|] and E731 ][Tj| experiments and recent results 



from KTeV |14| and NA48 [15]. The output parameters of the performed K — > 2ix analysis 
are the factors B\, and -B5 for a fixed value of B%. 

The dependence of Sj-factors on different theoretical uncertainties and experimental 
errors of various input parameters is investigated by applying the "Gaussian" method. 
Using Wilson coefficients derived in |J in various regularization schemes (LO, NDR, HV) 
for different values of the QCD scale A^L. ; we calculated the probability density distributions 
for Sj-factors obtained by using Gaussian distribution for all input parameters with their 
errors. As an example, the probability densities for the parameters B\, B±, B 5 calculated 
with B s = 1 and A^L = 325 MeV are shown in Figure [1|. We have shown the necessity 
for a rather large gluonic penguin contribution to describe the recently confirmed large 
experimental e value (the factor B§ is found well above 1). Figure |2| shows the correlations 
between B5 and Bg calculated for central values of all input parameters. From this figure 
one can see that even for B 8 = values of B 5 > 2 are necessary to explain the large value 
of e /e. It should be emphasized, that for even larger values of B 5 , the contributions of 
nonpenguin operators to the AI = 1/2 amplitude are still dominating (see Figure []). The 
large B\ and B§ values may be a hint that the long-distance contributions, especially to 
AI = 1/2 amplitudes, are still not completely understood. An analogous conclusion has 
been drawn in || , where possible effects from physics beyond the Standard Model are also 
discussed. 

Finally, we present our predictions for the CP-asymmetry of linear slope parameters in 
the — ► 37r Dalitz plot. These predictions are based on a new calculation of K — > 37r 
amplitudes at 0(p 6 ) within the same effective lagrangian approach. The obtained K —>■ 3tt 
amplitudes include the same theoretical uncertainties as in case of the K — » 2n analysis. 
The values of Bi, B4, B 5 fixed from the K — » 2n analysis are used as phenomenological 
input to the K — > 3n estimates which have been performed in a self-consistent way by the 
Gaussian method. 

The linear slope parameter g of the Dalitz plot for K — > 3tt decays is defined through 
the expansion of the decay probability over the kinematic variables 

\T(K -> 3tt)| 2 oc 1+gY + ... 

where Y is a Dalitz variable: 

v S3 - sp 2 m\ 2 

mi 3 
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Figure 2: Correlations between parameters B 5 and B 8 = 1. 
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Figure 3: Probability density distributions for the relative contribution of penguin operators 
to the AI = 1/2 amplitude. 
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Figure 4: Probability density distributions for the CP-asymmetry of linear slope parameters 
of K ± — > 7r ± 7r ± 7r T and K ± — > n 7r TT ± decays. 
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and the notation belongs to the non-equivalent pion in the decays K — > 7r ± 7r ± 7r =F 
and K — > 7r°7r 7r ± . Direct CP-violation leads to a charge asymmetry of the linear slope 
parameter, 



± _ _,_ g (K+ ^ 3n) - g(K- ^ 3n) 



Ag(K -> 3tt 



g(K+ -> 3tt) + flr(if- -> 3tt) 



In Figure |4] we show the probability density distributions for — > 7r ± 7T =t: 7r T and ^FtT^ — > 
Tr ^ ^ 111 decays calculated with Bg — 1 and = 325 MeV. Upper and low bounds for 

Ag + -i and Agoo+ for different values of A^g^ in LO, NDR and HV regularization schemes 

(B$ = 1) obtained by the Gaussian method are shown in table 1. The limits without 
brackets correspond to the confidence level of 68% while the limits in brackets - to the 
confidence level of 95%. Summarizing these results, we have obtained the following upper 
and lower bounds for the charge symmetries of the linear slope parameter: 

2.1 < A# ++ _-10 4 < 6.2, 1.4 < A# 00+ -10 4 < 3.5 with CL=68%; 

1.4 < Ag ++ _ • 10 4 < 10.4, 1.0 < Ag 00+ • 10 4 < 5.1 with CL=95%. 



Table 1. Upper and low bounds for Ag ++ _ and Ag 00+ (in units 10" 
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HV 
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215 


2.5 


5.9 


2.1 


4.7 


2.6 


6.0 






(1.8 


10.1) 


(1.5 


7.3 ) 


( 1.9 


10.2 ) 




325 


2.6 


6.0 


2.2 


4.7 


2.6 


5.8 






(1.9 


10.2) 


(1.5 


7.2 ) 


(1.8 


9.7 ) 




435 


2.7 


6.2 


2.2 


4.7 


2.6 


5.8 






( 1.9 


10.4) 


(1.4 


7.0) 


(1.8 


9.5 ) 




215 


1.5 


3.3 


1.4 


2.9 


1.6 


3.4 






(1.1 


4.8 ) 


( 1.0 


4.2 ) 


(1.2 


5.1) 


A#oo+ 


325 


1.6 


3.4 


1.4 


3.0 


1.6 


3.3 






(1.2 


5.0 ) 


(1.0 
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(1.2 
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1.7 


3.5 


1.5 


3.1 


1.7 


3.4 






(1.2 


5.1) 


( 1.0 


4.3 ) 


(1.2 


5.0 ) 



With some experimental updates and theoretical refinements our new estimates confirm 
the dynamical enhancement mechanism for the charge asymmetry Ag by higher order 
contributions in the effective chiral lagrangian approach which was first observed in ||. 
The predicted slope parameter asymmetry, although small, may be in reach of current high 



statistics experiments JT6[. Already with lower statistics, new measurements of quadratic 
slope parameters of K — » 3n decays, including the neutral channels, would lead to the 
improved theoretical understanding of the nonperturbative part of nonleptonic kaon decay 
dynamics. 
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